Effect of pulse error accumulation on dynamical decoupling of the electron spins of 

phosphorus donors in silicon 
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Dynamical decoupling (DD) is an efficient tool for preserving quantum coherence in solid-state spin 
systems. However, the imperfections of real pulses can ruin the performance of long DD sequences. 
We investigate the accumulation and compensation of different pulse errors in DD using the electron 
spins of phosphorus donors in silicon as a test system. We study periodic DD sequences (PDD) based 
on spin rotations about two perpendicular axes, and their concatenated and symmetrized versions. 
We show that pulse errors may quickly destroy some spin states, but maintain other states with high 
fidelity over long times. Pulse sequences based on spin rotations about x and y axes outperform those 
based on x and z axes due to the accumulation of pulse errors. Concatenation provides an efficient 
way to suppress the impact of pulse errors, and can maintain high fidelity for all spin components: 
pulse errors do not accumulate (to first order) as the concatenation level increases, despite the 
exponential increase in the number of pulses. Our theoretical model gives a clear qualitative picture 
of the error accumulation, and produces results in quantitative agreement with the experiments. 

PACS numbers: 03.67.Pp, 03.65.Yz, 76.30.-v 
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I. INTRODUCTION 

The state of a closed quantum system can be repre- 
sented as a coherent superposition of the basis states, 
where the phases and the amplitudes of the superposi- 
tion coefficients specify the state. Coherent superposi- 
tions of two-level quantum systems (qubits) can be clev- 
erly employed for high resolution nuclear magnetic res 
onance (NMR) and electronic spin resonance (ESR)^ 
highly sensitive magnetometers^— quantum information 
processing and quantum computation* 7 ^ However, any 
real quantum system interacts with its environment. The 
interaction leads to the decay of the coherence between 
the system's basis states, and the decohered qubit loses 
its useful properties. Different ways have been developed 
to mitigate decoherence£~— Among others, the dynam- 
ical decoupling (DD) approach looks extremely promis- 
ing. DD is based on the ideas that underlie the spin echo 
effect i 13 ' 14 A specially designed sequence of pulses is ap- 
plied to the system, and modifies the system-environment 
coupling in such a way that the impact of the environ- 
ment is averaged outj 15 ' 16 A large number of efficient 
DD sequences have been devised in the context of high- 
resolution NMR.— Experimental techniques for produc- 
ing the pulses are well developed, especially in the ar- 
eas of NMR and ESR. Furthermore, the threshold re- 
quirements for application of DD are very modest. All 
these advances make DD an appealing strategy in de- 
coherence suppression, and DD is actively studied now, 
both theoreticall y 16 ' 18 " — and experimentally^ 6 - - — 

Among different DD schemes, in this paper we focus on 
DD based on periodic structure, which has been broadly 
studied and implemented. A basic DD scheme is the pe- 



riodic dynamical decoupling (PDD), in which pulses are 
repeatedly applied to the system with equal inter-pulse 
delay 15 ' 16 To suppress the effect of system-environment 
interaction more efficiently, symmetrized versions of peri- 
odic dynamical decoupling ) 17 ' 71 and the concatenated dy- 
namical decoupling (CDD) have been proposed] 18 ' 23 ' 24 ' 38 
Performance of DD schemes can be different depending 
on the noise spectrum of the bath, and has been exten- 
sively studied for different qubit systemsi 20 ' 30 ' 31 ' 34-36 ' 72 

Besides the bath dynamics, real qubits are subject also 
to errors in the control pulses, which are often system- 
atic, being caused by the instrumental imperfections. 
Many research efforts have been devoted to achieving 
single pulses with small errors, or alleviating the ef- 
fect of finite pulse duration, including composite pulses, 
soft pulses and Eulcrian firli Yet implemen- 

tations of DD in QIP involve a large number of pulses 
applied over long times, so that even small pulse imper- 
fections could accumulate and seriously affect the decou- 
pling fidelity. A comprehensive understanding of how 
the systematic pulse errors accumulate in different DD 
sequences is highly desirable. The effect of pulse er- 
rors induced by the finite pulse duration has been stud- 
ied in detai h 30 ' 38 ' 66 ' 68 The effect of systematic errors in 
the rotation axis and angle on CDD has been studied 
theoretically^ and experimentally ) 34 ' 67 and has also been 
discussed for dynamically corrected gates j 6 - 6 - Moreover, 
for many traditional NMR/ESR experiments, the state 
to preserve is known (e.g. prepared along a certain di- 
rection) and the decoupling sequence can be chosen to 
best maintain this state. The Carr-Purcell-Mciboom- 
Gilli%>ulse sequence, where 7r-pulses are applied to the 
system along the £-axis, preserves the spin state along 
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the cc-axis while the state along the y-axis is destroyed 
by the accumulation of pulse errors. Quantum informa- 
tion processing, however, requires any unknown state of 
the qubit to be preserved, so DD must preserve all com- 
ponents of the qubit state. 

In this paper we analyze, both theoretically and ex- 
perimentally, the performance of several DD protocols in 
the presence of pulse imperfections, for different initial 
states. As a testbed, we use a macroscopic ensemble of 
electron spins of 31 P donors in silicon, which constitute 
a promising system for studying many fundamental as- 
pects of QIPM~— Instead of errors in a single pulse, we 
focus on the accumulation and compensation of system- 
atic pulse errors on a scale of the whole DD sequence. 
We study the performance of two-axis PDD, in which 
the 7r-pulses are applied along two perpendicular axes al- 
ternately, as well as its symmetrized and concatenated 
versions. Decoupling fidelities for spin states along x, y 
and z axes are examined. 

We find that, in PDD, certain types of pulse errors 
accumulated during the first half period of the sequence 
are balanced out during the second half, hence have little 
affect on the decoupling fidelity. Sequences with pulse 
axes along x and y alternately (XY PDD) are mainly 
affected by the in-plane component of the errors in pulse 
axes, while sequences with pulses along x and z axes 
(XZ PDD) are sensitive to multiple types of pulse error. 
Such analysis could serve as a guide for choosing suitable 
decoupling sequences according to specific experimental 
situation. 

We also find that some spin states arc quickly de- 
stroyed in the decoupling process, while other compo- 
nents are maintained with high fidelities over long times. 
Such a preservation of a particular spin component does 
not imply good performance of the decoupling sequence, 
but oppositely, results from the accumulation of the pulse 
errors over long tiniest E.g., XY PDD preserves the spin 
component along the z-axis, while XZ PDD preserves the 
y component. XY-bascd DD sequences are found to out- 
perform the XZ-based sequences. We analyze the depen- 
dence of DD sequences on the accumulation of pulse er- 
rors, and explain all these effects both qualitatively, with 
a simplified analytical model and quantitatively, using 
numerical simulations. 

We also find that the concatenated DD protocols 
(CDD) with two-axis control exhibit excellent immunity 
to pulse errors, and preserve all spin components. The 
robustness of CDD against systematic pulse errors has 
been shown theoretically earlier^ Our results provide 
experimental confirmation of this feature of CDD: in the 
evolution operator the pulse errors do not accumulate 
to first order, in spite of the exponential growth of the 
number of pulses. 

The rest of the paper is organized as follows. In Sec.lTIl 
we introduce the system and the DD sequences to be 
studied. We present the theoretical model, and describe 
the experiments. In Sec. lIIIi we present the results of an- 
alytical studies and numerical simulations, and compare 



them with the experimental data. Conclusions are given 
in Sec. \N\ 



II. THE SYSTEM AND THE DD PROTOCOLS 
A. The phosphorus doped silicon system 

Electron spins of P donors in silicon show long 
relaxation time s 44 ' 45 and coherence times,— and 
therefore constitute promising candidates for QIP 
application s 39 ' 40 ' 47 ' 48 and for studying the basic prob- 
lems of DDi 26 ' 40 Advanced ESR techniques can be used 
to manipulate the spin states of the P donorsi 46 ' 49 In 
our experiments, isotopically purified bulk silicon sam- 
ples were used; 5 - - with 29 Si concentration ~ 800 ppm, 
and the doping density of phosphorus ~ 5 x 10 14 cm -3 . 
Experiments were performed at a static (quantizing) field 
of 3500 Gauss, and a temperature of 8 K. The longitu- 
dinal relaxation time of the electron spin of P is T\ = 
25 ms limited by the two-phonon Orbach processi 46 ' 51 
The transverse relaxation time determined by spin echo 
experiments is T2 = 4.6 ms, limited by instantaneous 
diffusion! 46 ' 52 ' 53 The "true" T2 of isolated donors extrap- 
olates to about 60 ms, e.g. after suppressing the instan- 
taneous diffusion^ 6 - 

Dynamical decoupling works well when the inter-pulse 
delay between the pulses is short compared to the typical 
time scale of the noise of the bath. In our experiments, 
the inter-pulse delay time is r = 11 /us. We now examine 
the time scales of possible interactions with the P electron 
spins. 

At donor density ~ 5 x 10 14 cm -3 , the typical mutual 
flip-flop time of the P electron spins due to the dipolar 
interaction is of the order of 100 ms, which is much longer 
than the inter-pulse delay in the DD and the total dura- 
tion of the experiment. We therefore can regard the P 
electron spins as independent of each other. 

The hypcrfinc coupling between the donor's electron 
spin and the donor's 31 P nuclear spin is rather large 
(~ 100 MHz)j^ and only one hyperfine line is excited in 
our experiments (i.e. we are working only within the sub- 
space with a fixed z-projection of the 31 P nuclear spin). 
Thus, the presence of the 31 P nuclear spin only slightly 
renormalizes the resonance frequency of the donor's elec- 
tron spin, and can be neglected. 

Spectral diffusion induced by the nuclear spins is a 
major source of decoherence for Si:P system and other 
dopants in semiconductors i 5 - 5 -"— In Si:P system, the P 
donors are coupled to their surrounding 29 Si nuclear spins 
mainly via contact hyperfine interaction. The anisotropic 
corrections caused by the admixture of the p-states and 
the dipolar contribution to the hyperfine coupling, are 
smalli 5 ^ The dynamics of the 29 Si nuclear spin bath, ei- 
ther due to the flip-flops induced by the intra-bath dipo- 
lar interaction or due to the electron-mediated virtual 
spin flips, could decohere the P spins ^ r 58 i 60 However, 
in our experiment the concentration of 29 Si is low (^800 
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ppm) and the quantizing field Bq = 3500 Gauss, which 
we take as directed along the z-axis, is huge compared to 
the hyperfine coupling energy scale. Both types of flip- 
flops of the nuclear spins are thus greatly suppressed on 
the time scale of r and on the timescale of the total ex- 
periment. The characteristic spectral diffusion time has 
been measured to be about 20 ms at 800 ppm of 29 Si, fili 
which is larger than the inter-pulse delay by four orders 
of magnitude. Therefore the 29 Si nuclear spins can be 
treated as static, contributing only an extra static field 
acting upon the P spins along the z-axis. Moreover, since 
the separation between the donor centers is very large, 
the 29 Si nuclei which are sufficiently strongly coupled to 
one center interact very weakly with other donor centers 
and each P electron spin can be considered as coupled to 
its own nuclear spin bath. 

The quantizing magnetic field of the magnet also fluc- 
tuates. The fluctuations have a broad spectrum but the 
noise power quickly decays as l// 2 to become negligible 
at frequencies higher than 1 kHz. The noise correlation 
time is therefore much larger than r and we can treat 
the quantizing field as static but inhomogeneous over the 
sample. 

Therefore, the system under study is an ensemble of 
independent electron spins. Each spin feels a different 
static quantizing field, and an extra static field con- 
tributed by the hyperfine interaction with the 29 Si nu- 
clear spins. The Hamiltonian describing a single phos- 
phorus electron spin S is then 



H = ui e S z + j e B n S z 



(1) 



where w e is the Zceman frequency of the electron spin in 
the static quantizing field, 7 e is the gyromagnetic ratio 
for the P electron, and B n is the overall effective field 
acting on S due to interaction with the bath. Here and 
below we adopt the system of units with h = 1. Tak- 
ing into account the spatial inhomogeneity of the static 
quantizing field, we represent cu e = cu e Q + 5uj e , where w e o 
is the average over the sample, i.e. the frequency corre- 
sponding to the center of the ESR line. By performing 
the standard rotating-framc transformation ; 13 ' 43 uj e o is 
eliminated and the Hamiltonian (JT|) is transformed into 



H = le BS\ 



(2) 



where B = 5uj e /"/ e + B n is the total field acting upon S in 
the rotating frame. Accordingly, all analysis below is per- 
formed in the rotating frame. The distribution function 
of this field is determined in part by the initial density 
matrix of the bath, and in part by the inhomogeneity of 
u! e over the sample. The statistical properties of the field 
B determine the ESR lincshape. Our measurements of 
the free induction decay show that the ESR line has a 
Gaussian shape, hence 



P(B) = 



1 



:ex P [-B 2 /(2& 2 )] 



(3) 



with variance b — 50 mG, which includes both the rms 
coupling between the spins of the P donors and the bath, 



and the inhomogeneity of the external quantizing field. 
This experimental fact will be the starting point of the 
further analysis. 

To sum up, on the time scale of the inter-pulse delay 
of DD (tens of microseconds), the electron spins effec- 
tively live in a static inhomogeneous background mag- 
netic field. The dephasing time T 2 * ~ 1.6 /xs of the P 
spins due to such inhomogeneity is much smaller than 
Ti . Sequences of microwave- frequency pulses are applied 
to decouple the P electron spins from the overall back- 
ground magnetic field noise, and refocus the dephasing. 
The electron spin state is expected to be preserved by 
DD against the dephasing on a time scale of the order of 
T 2 , where instantaneous diffusion becomes important. 



B. Experimental setup 

Our experiments were performed on a Brukcr Elexsys 
580 spectrometer with specially modified software to al- 
low for generating large numbers (over 1000 in some ex- 
periments) of microwave pulses, and a 20 watt contin- 
uous wave solid state microwave power amplifier (Am- 
plifier Research) which maintained microwave phase sta- 
bility over the long pulse sequences. Experiments were 
performed at an X-band microwave frequency of 9.8 GHz 
and a static (quantizing) magnetic field of 3500 Gauss, at 
temperature 8 K. The typical duration of a 7r-pulse was 
0.18 /xs. In most of the experiments the delay between 
refocusing pulses was set to r = 11 /xs, short compared to 
the characteristic times of all sources of noise known in 
our system (see the detailed discussion in Section III Al) . 

For the donors in our Si sample the measured echo sig- 
nal decays showed an excessive "phase noise" developing 
at longer times (for example, at t > 1 ms in a Hahn 
echo experiment) arising from the fluctuating magnetic 
field Bq of the magnet.— The power spectrum of the field 
noise (Gauss 2 • Hz -1 ) varies approximately as l// 2 with 
an amplitude ~ 0.5 mG-Hz -1 ' 2 at 10 Hz. In the experi- 
ments we used a large Si sample, with enough donors to 
enable the acquisition of an echo without signal averag- 
ing. Thus the field noise could be eliminated with mag- 
nitude detection, e.g. squaring and adding the in-phase 
and quadrature components of the echo signal.— This 
magnitude detection approach was used for an accurate 
extraction of the Ti relaxation times from the measured 
decays. 

The large Si sample also means that the mi- 
crowave magnetic field (B p ) is not homogeneous over 
its entire volume in the pulsed ESR resonator (a 
Brukcr dielectrically-loadcd cylindrical cavity 4118X- 
MD5), leading to a systematic rotation angle error. The 
extent of B p inhomogeneity and the resulting rotation er- 
rors in our spectrometer setup have been characterized in 
our previous publications.— ^ The B p variation over the 
sample volume can be as large as 10-20% depending on 
the sample dimensions. On the other hand, the rotation 
axis errors (e.g. relative phase errors between 7Tx and 
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7Ty pulses) can be reduced to a sub-degree level using 
a calibration technique based on phase error amplifica- 
tion (SPAM),— in which the phase errors accumulate in 
a way which is similar to how the flip angle errors accu- 
mulate in the CP sequences, and hence the orthogonality 
of the rotation axes in two pulse-forming channels can be 
measured with high precision. 



C. Dynamical decoupling protocols 

In a spin echo experiment, the in-plane magnetiza- 
tion of a spin ensemble decays due to the precession of 
the spins under an inhomogeneous magnetic field, and 
this dephasing can be refocused by applying a 7r-pulsc 
to all the spins midway during the evolution. Pulse se- 
quences like Carr-Purcell-Meiboom-Gill (CPMG) refocus 
the dephased ensemble magnetization by applying a train 
of equally spaced pulses, each rotating the spin by it 
around the x-axis of the rotating coordinate frame. Be- 
low we denote such pulses as 7rx, and other pulses are 
denoted similarly (e.g. a (7r/2)Y-pulse is the pulse rotat- 
ing the spins by 7r/2 around the y-axis). However, the 
single-axis CPMG protocol has an important drawback. 
When the electron spins are prepared along the y-axis 
(by applying a preparatory (ir/2)x pulse), the small er- 
rors inevitably present in real DD pulses accumulate very 
quickly in the course of the decoupling experiment, and 
destroy the performance of the decoupling protocol, like 
in the original Carr-Purcell pulse sequence. This prob- 
lem is absent for CPMG in case of the initial states along 
x-axis, but the goal of DD is to preserve all components 
of the spin. Therefore, more general DD protocols must 
be considered. A large number of decoupling sequences 
have been developed, aiming at maintaining the state of 
the central spin S. Design and analysis of DD schemes 
are often based on the Magnus expansion (ME), which 
is a cumulant expansion of the evolution operator of the 
whole system (the qubit and the bath). DD intervenes 
the evolution and eliminates the unwanted system-bath 
couplings in expansion terms of the ME. Besides the DD 
sequences based on a periodic pulse structure like PDDi& 
and its concatenated version CDD^ aperiodic structure 
like UDD^ QDD^ etc. were developed recently. In this 
paper we focus on the decoupling sequences with periodic 
structure. 

One of the simplest protocols is the periodic dynamical 
decoupling (PDD),— where the central spin is subjected 
to a train of equidistant 7r-pulses applied along different 
axes. The unitary operators, which describe the con- 
trol pulses acting on the central spin, are taken from a 
group cyclically, starting from the identity element I. A 
typical sequence based on the group Q = {1 7 a x ,a v ,ct z } 
where a x , a y and a z are Pauli matrices has a period 
(I-d-I) (7r x -d-7Tx) (7r z -d-7r z ) (7rY-d-7r Y ) . Here d indi- 
cates the delay between pulses with time duration r. By 
virtue of the operator algebra of the Pauli matrices, the 



period can be reduced to 

d — 7T X — d — 7T Y — d — 7T X — d — 7T Y , (4) 

This two-axis (XY) based decoupling scheme is termed 
XY PDD. When the inter-pulse delay time r is short 
compared to the inverse of the cutoff frequency, w c , of the 
spectral density function of the bath, the resulting evolu- 
tion operator of the system is equivalent to the identity 
up to the first order in w c t, i.e. has a form 1 + C>(w^t 2 ). 
In this paper we also study an alternative two-axis de- 
coupling sequence, XZ PDD. Its period 

d — 7T X — d — 7T Z — d — 7T X — d — 7T Z . (5) 

can also be drawn from the group Q. For ideal pulses 
(e.g., when the bath is completely static and pulse errors 
are absent), XY- and XZ-based DD are equivalent to each 
other. 

SDD is a symmetrized version of PDD. The period of 
SDD sequence is twice as long as that of PDD and is 
symmetric with respect to the middle i 16 ' 17 ! 71 The period 
for a XY-based SDD is 

d — 7T X — d — 7Tf — d — Tlx — d — 7Ty 

— 7Ty — d — 7T X — d — 7Ty — d — 7T X — d . (6) 

Such symmctrization eliminates all even terms in the ME 
systematically, leaving the system-bath coupling in the 
third and higher odd order terms. 

A concatenated version of periodic decoupling, CDD, 
was developed to decouple the system from the bath to 
higher orders in ME^ CDD at level one is just PDD, 
and sequences for higher levels are built up by recursively 
nesting the lower level sequence within itself. E.g., for the 
XY-based sequence, the structure for n-th level CDD, 
CDD n , is 

CDD n _ t -7T X -CDD n _ t -7T Y -CDD n _ t -7T X -CDD n _ I 7Ty (7) 

In ME, CDD n eliminates the interaction between the sys- 
tem and the bath up to the n-th order. Note here the 
number of pulses in CDD increases approximately as 4™. 

As is customary in NMR/ESR experiments, the ini- 
tial state of the electron spin is pseudo-purei 13 i 43 The 
density matrix for the electron spin polarized along the z- 
axis in the high-temperature approximation has the form 
po = |l + C,S Z where C is small. Since the identity ma- 
trix 1 is not affected by unitary evolution and gives no 
contribution to the signal, we can regard the electron 
spin as being in a pure state with S z = 1/2 (pseudo-pure 
state) i 13 ' 43 The in- plane initial states of the electron spin 
prepared by (tt/2) -pulses at the beginning of the experi- 
ment then can also be considered as pure. 

For a single electron spin S, we characterize the per- 
formance of a given DD protocol using the survival prob- 
ability of the state (input-output fidelity), which in our 
case can be written in the form Tr[p(0)p(t)], where p(0) 
is the initial state of the electron spin and p(t) is the 
reduced density matrix of the electron spin at the end 
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of DD. By virtue of the relation p(t) = 1/2 + (S x )a x + 
(S y )a y + (S z )a z , the fidelity is equivalent to 

F%(t) = 2Tr{S a (t)S a (0)} = 2Tr[p a (t)S a ] , (8) 

and is characterized by the average spin projections 
2{S X ), 2(S V ), and 2(S Z ). Here a = x,y,z denotes the 
three initial states of the electron spin oriented along the 
axes x, y and z respectively, and p a (t) is the correspond- 
ing density matrix of the spin at time t. The fidelity char- 
acterizes how well the initial state of the electron spin is 
protected by the DD sequence. F^(t) is then averaged 
over the ensemble as 

F a (t) = {F%(t)), (9) 

where (•) denotes the ensemble average. The quantity 
F a (t) is the measure of the performance of DD in our 
study. Of course, this fidelity may, and in fact does, 
strongly depend on the specific initial state. 

D. Analysis and model of pulse errors 

In ESR experiments pulse errors can be greatly re- 
duced, e.g. by using composite pulses^ Still, in many 
experiments pulse errors remain an important (and some- 
times a major) factor which limit the performance of the 
decoupling sequence. We now undertake a detailed anal- 
ysis of the possible pulse errors in our experiment, and 
arrive at a model to account for the error effect. 

Since the duration of a 7r-pulse, 0.18 ps, is small com- 
pared to the inverse ESR linewidth (~ T 2 *), we treat the 
pulses as instantaneous unitary rotations. The validity 
of this approximation has been confirmed numerically, by 
modeling the influence of the dephasing field B (Eq. fl2J)) 
during the pulses. Taking into account the errors in both 
the rotation axis and the rotation angle, the operators 
for the 7r-pulses have forms 

U x = expH(7r + e x )(S-n)] (10) 
Uy = exp + e y )(S ■ m)] , 

where n = ( yl 1 — n 2 — n 2 , n y , n z ) is the actual ro- 
tation axis for a nominal 7rx-pulse, and m = 
(m x , yl — rn 2 — m 2 , , m z ) is the actual rotation axis for 
a nominal 7TY-pulsc. Small parameters n y , n z , m x , m z 
characterize the error in the rotation axes, and e x and e y 
characterize the error in the rotation angles. 

1. Error sources 

An ideal 7r-pulse generated by a pulse field in exact 
resonance with the Larmor frequency of the central spin 
realizes a rotation of the central spin by an angle ir. In 
the rotating frame, the pulse field B p has a step-like form 
in time, i.e., the value of the field is zero before and after 



the pulse, and constant during the pulse. The duration t p 
of the pulse is determined by t p ^ e B p — ix. In this study, 
we consider the following four sources of pulse error. 

Due to the inhomogeneity in the static magnetic field 
over the sample, the pulse field is not in exact resonance 
for every electron spin. In the coordinate frame rotat- 
ing at angular frequency w e o, the detuning, B, manifests 
itself as a non-zero magnetic field along the z-axis. As 
a result, the rotation axis deviates from the intended x 

or y directions, and the rotation angle i p 7 e yj B 2 , + B 2 is 

different from tt. 

In our experiment, the coherence signal was obtained 
by a single-shot measurement from a macroscopically 
large Si sample. The byproduct of a large sample is that 
B p is not homogeneous over the entire sample volume. 
To model this effect in a simple but physically meaning- 
ful way, we assume a one-dimensional model, where B p 
varies only along some axis I (e.g. along the resonator 
axis), and the sample position is optimized near the maxi- 
mum of the driving field, where B p depends quadratically 
on I: 

B p {l) = B p + AB p [l-3l 2 /d 2 ], (11) 

where B p is the average of B p (l) over the sample, and 
AB p quantifies the amplitude of the error. The origin 
of the /-axis is at the center of the sample, and 2d is the 
sample length along the Z-axis. In the experiment, B p can 
be well tuned such that on average an exact 7r-pulse can 
be expected, i.e. t p j e B p = tt. The error in the rotation 
angle ir + e(l) for the spin located at position I then arises 
from the AB p term 

e(0 = t ple AB p (l-3l 2 /d 2 ) . (12) 

The corresponding distribution function for the rotation 
angle error is 

P(6) = (l/2 eo )[3(l-£/£o)r 1/2 . (13) 

and — 2eo < e < eo We further assume the angle errors 
for 7Tx- and Try- pulses to be the same e x = e y = e. 

The imperfectly rectangular shape of the actual pulses 
also introduces pulse errors. These errors mainly come 
from the leading and trailing edges of the pulse. In ESR 
experiments, the edges may constitute about 10% of the 
total pulse time, and could have a noticeable influence. 
At the pulse edges, the amplitude and the phase of the 
driving field arc ill-defined, introducing errors to both the 
rotation axis and the rotation angle. Furthermore, the 
magnitude of these transient errors depends on the offset 
frequency, B, of the spin and therefore the pulse errors 
have a non-uniform distribution over the sample. To ac- 
count for the effect of these complex errors, we introduce 
the rotation axis errors n y , m x and n z , m z . Similar to the 
arguments used in Eq. (|13p . the axis error, n z , is drawn 
from the probability distributions 

P{n z ) = (l/2n )[3(l-njn )r 1/2 , (14) 



6 



where — 2uq < n z < no, and constant no quantifies the 
amplitude of the error. Similarly, we assume the axis er- 
rors of 7TY-pulses have the same probability distribution 
as Eq. (|14p . and assume m z = n z for all spins. The value 
for the in-plane component of the rotation axis errors, 
m x and n y , are extremely small in experiments as dis- 
cussed in Section Hi B[ so we take them to be zero in the 
simulations, but keep the symbols for the analysis. 

Another possible pulse error is the imperfection in the 
relative phases of the 7rx and 7TY-pulses. In the presence 
of such errors, if we assume that the axis for 7rx-pulse 
is perfectly aligned along the x-axis, the actual rotation 
axis of the 7rY-pulse may have a non-zero x-component. 
However, by using standard phase calibration techniques, 
this in-plane axis error can be reduced to a sub-degree 
level and is negligible in our experiments. 

Summarizing, we treat the pulses as instantaneous ro- 
tations described by Eq. (|10|) . The pulse errors e x = e y = 
e and n z = m z are drawn from distributions Eq. (|13p 
and Eq. fH|) . respectively. Note that in our model for 
pulse errors, eo and no are the only adjustable param- 
eters which are determined from experiments. With a 
fixed set of values for eo and no, good simulation results 
are achieved in quantitative agreement with the experi- 
ments for all seven different DD protocols (XY and XZ 
PDD and CDD. Symmetrized DD sequences were not 
performed experimentally) . 

2. Simulation details 

Time evolution of the electron spin S is simulated for 
initial states along axes x, y and z. We label these three 
initial states respectively as \ip x ), \ipy) an( A \ipz), and the 
fidelities as F x , F y and F z . Simulations are performed 
in the rotating frame with the frequency at the center of 
the ESR line, u; e0 . For each single spin, the pulses are 
implemented as described by Eq. ([TO]) , and the evolution 
operator for the inter-pulse delay is Uq = exp [— iHt] with 
the Hamiltonian H given by Eq. @. The static field B 
felt by spin S is drawn from the Gaussian distribution, 
Eq. ([3]). After application of the decoupling sequence, 
the input-output fidelities, F^(t), are calculated. The 
fidelities are averaged over ~ 10 4 realizations. The time 
delay between pulses is constant r = 11 /is for all DD 
sequences as in the experiments. The values of the error 
parameters in Eqs. (|13[) . (fH)) arc taken as eo = 0.3 (7.5°) 
and no = —0.12 (—3.5°), which provides results that fit 
well with the experiments. 



III. RESULTS AND DISCUSSION 

Experimental and numerical results for different DD 
protocols are shown in Figure [1] through Figure [4] The 
simulation results agree well with the experimental data, 
indicating that the error model captures the essential fea- 
tures of the system. We present the results for decoupling 



sequences based on XY and XZ pulses for comparison. 
In the experiment, a 7rz pulse was implemented by a pair 
of subsequent rotations about the x and y axes, 7rx7TY- 
This substitution was carried out in experiments for all 
the XZ based sequences. Since the bath of the electron 
spins is static in our study, for ideal pulses (without pulse 
errors), XY- and XZ- based sequences are exactly equiv- 
alent. However, as we show below, in the presence of 
pulse errors they behave very differently in preserving 
the quantum states of the electron spin . 

A. PDD 

1. XY-based sequence 

Figure Q] (a) shows the fidelities of the electron spin 
states as functions of the number of periods (cycles) for 
XY PDD. Fidelities for different initial states exhibit 
different decay behaviors. F x and F y exhibit apparent 
decays. F z , on the contrary, decays very little. Simi- 
lar phenomenon in DD have been discussed in Ref. [23| . 
where one state of the system is less sensitive to the de- 
coherence induced by the spin-bath interactions than the 
other states. In our study, in the presence of pulse im- 
perfections, particular initial states of the electron are 
preserved, i.e., survive for a long time in the DD process, 
while other states are quickly destroyed by the accumu- 
lation of pulse errors. 

The electron spins of P are independent of each other, 
and are coupled only to external magnetic fields (static 
field along the z-axis, and a time-varying pulse fields 
along all three axes). In this case, the time evolution of 
each spin is unitary, and can be described as a rotation: 

U = exp [-i0(S • a)] , (15) 

where a is the effective rotation axis and 9 is the rotation 
angle. In the following we will examine such rotations, 
viewed stroboscopically after each period of the DD se- 
quence. 

Consider a period of the XY PDD sequence Eq. ([?]), 
substituting the expressions for pulse errors, Eq. (|10[) . 
into the evolution operator 

^xypdd = UyU U x U U y U U x U , (16) 

and keeping only the zeroth and first order terms in the 
small parameters (e x , e y , m x , n y , m z , n z ), [/ XYPDD can 
be expressed in form of Eq. (fTS"]) with 

a = (0,0,-1) 

6 = 2n + S9 (17) 

where 

89 = 4(m x + n y ) . (18) 

That is, after one period, the spin is rotated about an 
axis close to — z by an angle 2ir + 89. For ideal pulses, 
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jjxypdd _ the rotation angle is = 2tt and the spin 
returns exactly to its initial state. The error in the rota- 
tion angle, SO, can be viewed as a dcphasing error in the 
plane perpendicular to a. After N periods, to first order 
in the pulse errors, the spin rotates about the axis a by an 
angle NSO. Since a is (approximately) along the z-axis, 
the spin operator S z (approximately) commutes with the 
evolution operator, Eq. ([T5J) . Therefore the initial state 
\ip z ) is an (approximate) eigenstate of the evolution op- 
erator, and thus does not evolve with the DD sequence, 
and is of course insensitive to pulse errors (pulse errors in 
high order could contribute though). On the other hand, 
the other two states \tp x ) and \ip y ) are not the eigenstates 
and for these states the error in the rotation angle SO will 
accumulate during the pulse sequence. 

Now we consider the electron spin ensemble. For 
each spin the time evolution is governed by a rotation, 
Eq. (|15j) , with rotation axis and angle given by Eq. (|17jl . 
Each spin sees different pulse errors, and the resultant 
effective rotation axis, a, and dephasing errors, SO, are 
slightly different from each other (for first or higher or- 
der pulse errors). As shown in Eq. (|17j) . the rotation 
axes for all spins are along z up to first order in the pulse 
errors. The initial state \ip z ) is then preserved for all 
spins in the ensemble, despite their different pulse errors. 
The ensemble-averaged fidelity, F z , hence remains high 
until spoiled by the accumulation of higher order errors. 
On the other hand, for states \ip x ) and \tp y ), the spins 
are initially in the plane perpendicular to the a (z)-axis. 
Different spins acquire different phase errors during the 
evolution in one period of DD. As a number of periods N 
increases, this error accumulates. After certain number 
of periods, the spin components spread out evenly in the 
x-y plane, and the ensemble-averaged fidelity decays to 
almost zero. 

Note here if we analyze the rotation operator after half 
a period d — 7r x — d — 7r Y , (which is the repeating unit 
of the pulse sequence), and express the corresponding 
evolution operator L/*™ D up to first order in the pulse 
errors, the rotation angle is 0' = n + 2(m x + n y ) and the 
axis of the rotation, a', is 

a x = - — + n z cos (pd — — sin q>d 



m z - — cos <p A - n z sin <p d 



= -1 



(19) 



where (f>d = "f e Br is the phase accumulated during the 
inter-pulse delay. It is noteworthy that unlike the rota- 
tion axis for a full period in Eq. (fTT)) , the rotation axis a' 
involves pulse errors to first order. That is, some pulse 
errors accumulated during the first half period are bal- 
anced out during the second half. Since a full period is a 
repetition of two half periods, we certainly have 



U 1 



= « D XaT f DD = cx P H2^S-5')], 



(20) 



and the rotation axis is still the same and the rota- 
tion angle is twice larger. Although the expressions 



for a and a' are different, there is no discrepancy here. 
Since exp [-i20'(S • a')] = cos6»' - i/2sin0'(S • a') and 
sin6>' = — 2(m x + n y ), first order pulse errors in a' be- 
come higher order in (/ XYPDD . Therefore, up to fist order, 
L/xypdd can also be expressed as Eq. (|T7| with the effective 
axis, a, free of first order pulse errors. 

We next undertake a more quantitative examination of 
how the decays and the preservation happen to different 
spin components. For a single spin, S, the fidelity for an 
initial state, \ip a ) i a = x i 2/? z )i a ft cr N periods is 



2Tr [U- 



rPDD ps(0)(C/ XYPDD ) t ^° 
a 2 a ) cos (N 50) . 



(21) 



That is, the fidelity, if viewed stroboscopically (at the end 
of each DD period), exhibits oscillations at a frequency 
determined by the phase error, SO. For spin ensemble, 
after averaging over B and over the pulse errors (although 
the magnetic field B does not show up in Eq. (fT5|) . it 
contributes in higher order) , the cosine factor contributes 
to the decay, and the fidelity at long times saturates at 
a value independent of time (N) 



F a (N) 



(22) 



where (■) denotes the ensemble average. For \tp x ) and 
\ijjy), (a^) is small, hence the fidelities decay almost to 
zero. For \ip z ), since a z is almost one, the fidelity is 
preserved with its value close to 1. 

Such preservation of a particular spin state demon- 
strates a remarkable feature of the error accumulation in 
DD: a large fidelity for certain initial states does not guar- 
antee good performance of DD for other initial states, 
hence does not necessarily suggest a good DD protocol. 
This is similar to the error accumulation effect in the 
single- axis DD sequences. While \ip x ) is well preserved 
by the 7rx-pulsc sequence (CPMG), \ip y ) is quickly de- 
stroyed by the same sequence due to the accumulation of 
errors in the rotation angle (CP). Sensitivity of Uhrig 
DDi2 and CPMG to initial input state in the presence of 
pulse errors also has recently been studied experimentally 
in Rcf. 67j. Therefore, in experiments implementing 



F y and 



DD, the preservation of all three components F, 
F z should be checked, or full tomography should be per 
formed. 



2. XZ-based sequence 

The pulse sequence of XZ PDD preserves the spin com- 
ponent \tp y ), see Figure [T](b). Compared to the decaying 
fidelities {F x and F y ) in XY PDD, the fidelities F x and 
F z in XZ PDD decay much faster, with values below 0.2 
after barely 3 cycles (12 pulses). 

An analysis similar to that in IIII A II shows that the 
evolution operator for one period of XZ PDD {/ XZPDD = 
exp [— i(2ir + SO) (S ■ a)], up to first order in the pulse er- 
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FIG. 1: (color online). Fidelity as a function of the number of 
cycles, N, of PDD. (a) XY PDD, (b) XZ PDD for initial states 
\ipx), l^y), an d \i>z), respectively. Points are experimental 
results and lines are simulation results. 

rors, corresponds to 

a = (0,-1,0) 

89 = 2[— e y + e x sin^d + 2n z (l — cos^d)] ■ (23) 

The effective rotation axis a is close to —y, hence the 
state \ip v ) is preserved. 

Comparing Eq. to Eq. (|T7|). one can sec that XZ 
PDD involves rotation angle errors e x , e y and axis error 
n z to first order, while the XY-based sequence involves 
only m x and n y . Since in experiments, as we mentioned 
in Section TlI D[ the in-plane rotation components, m x and 
ii y , are extremely small and can be neglected, the devi- 
ation of the evolution operator of a period of XY PDD 
from identity is actually second order in the pulse errors 
while that of XZ PDD is first order. Furthermore, note 
that Eq. (f!?3"|) for XZ PDD is dependent on the magnetic 
field (in <j> d ) while Eq. Q2|) for XY PDD is not. When 
performing the ensemble average in Eq. (f!2~Tj) . which holds 
also for XZ PDD and in general for any periodic decou- 
pling sequence, the first order pulse error in 89 and the 
independence leads to a much faster decay of the cosine 
term than the second order pulse error does for XY PDD. 
This explains why XY PDD outperforms XZ PDD. 

Note that such inferior performance of XZ PDD com- 
pared to XY PDD does not arise from the replacement of 



7Tz by 7Tx7Ty in our experiment. The XZ-bascd sequence 
still performs worse even if 7Tz-pulses are applied directly. 
Taking the form of the evolution operator for 7Tz to be 
similar to Eq. (jTOJ) , with rotation axis p and rotation 
angle n + e z , Eq. (|23)l becomes 

a = (0,-1,0) 

89 = 2[-2p x + e x sin0 d - 2n z cos0 d ] , (24) 

where p x is the x-componcnt of the rotation axis for 
the 7Tz-pulsc. (Note, in the limit 4>d 0, we have 
89 = — i(p x + n z ), which is exactly the symmetric form of 
their counterpart for the XY-based sequence Eq. (fT7|) , as 
expected from the rotational symmetry of the system.) 
The evolution operator for XZ PDD hence contains mul- 
tiple first order pulse errors, similar to the case where 7Tz 
is replaced by itxtty- Therefore, given the experimental 
facts that the in-plane component of the pulse axis error, 
m x and n y , can be made extremely small while the other 
errors cannot be neglected, we can expect XZ PDD to 
perform worse than XY PDD. 

In summary, pulse errors play an important role in 
periodic DD. Each PDD sequence has an effective axis 
in the presence of pulse errors. During the decou- 
pling, viewed stroboscopically, the spins perform rota- 
tions about this axis. Error accumulation could destroy 
the initial state components perpendicular to the effec- 
tive rotation axis, while the component parallel to the 
axis is preserved. For negligible pulse axis errors, m x 
and n y , XY based decoupling sequence shows a better 
resistance to the pulse imperfections than the XZ-bascd 
sequence. 



B. CDD 

Next we examine concatenated sequences based on 
PDD. We measured the fidelity for a single cycle of XY 
and XZ based CDD, with concatenation levels 1 to 4. 
Correspondingly, the numbers of pulses are 4, 20, 84 and 
340, and the number of delays (all delays have the same 
duration r) are 4, 16, 64, and 256, respectively. Fidelities 
are evaluated after the application of the whole sequence. 

In XY CDD, as shown in Figure [5] (a), all three ini- 
tial states are preserved with excellent fidelities (> 99%). 
Numerical simulations based on our error model are in 
agreement with the experimental results. Note here, at 
level 4, the electron spin has been subjected to 340 pulses, 
and its quantum state is perfectly preserved with fidelity 
close to unity for 2.8 ms, which is of the order of T-2,. For 
many pulses and over a long time, no visible sign of the 
error accumulation appears. 

To better understand the excellent performance of XY 
CDD, we analyze error accumulation in the concatena- 
tion. The evolution operator for a period of the level-2 
CDD [/ XYCDD2 j U p to first order in the pulse errors, corre- 



9 




0.5 







1 1 1 1 

■i A 


NN. ▲ 


^ 


© 


(b) 





12 3 
Concatenation level 



FIG. 2: (color online). Concatenated DD. Fidelity as a func- 
tion of the concatenation level is shown for initial states \tpx), 
\ip y ) and \ipz)- Points are experimental results and lines are 
simulation results, (a) XY CDD. All three lines for states 
IV*}, \ip y ), and \ij} z ) overlap, (b) XZ CDD. 



For the first level concatenation (PDD) the evolution op- 
erator depends on the phase during the inter-pulse delay, 
4>d, i-e. depends on the factor 7 e i?r, but the higher lev- 
els are independent. This explains the dips at level- 1 
concatenation for F x and F z in Figure [2] (b) . Therefore, 
the concatenation eliminates the dependence of the error 
accumulation on the inter-pulse delay. 

The simulation results for XZ CDD are qualitatively 
in agreement with the experimental data, e.g. state \tp y ) 
is better preserved than the other two states, and there 
is a dip in the first-order concatenation for \ip x ). The 
quantitative difference may be due to other experimental 
details that we did not take into consideration in our 
error model. 



C. Symmetrized DD 

We also examined the symmetrized XY-based decou- 
pling sequence via numerical simulations and theoretical 
analysis. The symmetrized pulse sequence unit Eq. ((6]) 
is repeated periodically (XY-SDD). Figure [3] shows the 
simulation results. The fidelities of the spin states exhibit 
much slower overall decay compared to XY PDD. Indeed, 
thanks to the symmetrization, the first order pulse errors 
balance out in the evolution operator for a single period 
of SDD. Up to second order, the evolution operator reads 

(yXY-sDD = 1 + 2ie y (m x + n y )a x 

+2i(m x + n y ) (e x cos 4>d + 2n z sin </>d) a y . 



(28) 



sponds to a rotation with rotation axis and angle 
a = (0,0, -I) 

56 = 4(m x + n y ) . (25) 

This exactly repeats the rotation corresponding to XY 
PDD, Eq. (TIT)) . Further examining higher concatenation 
levels, we arrive at the following relations 

jjXYPDD _ jjXYCDD2 _ . . . _ ^jXYCDDn ^g) 

where [/ XYCDDl1 denotes the evolution operator for level-n 
concatenation. Up to first order in the pulse errors, the 
evolution operators are exactly the same for all concate- 
nation levels. Although the number of pulses increases 
exponentially with the concatenation level, the first or- 
der errors do not accumulate. Concatenation endows DD 
with the capability of self-correcting the pulse errors. 

XZ-based concatenated DD also outperforms its PDD 
counterpart, see Figure [5] (b). However, compared to XY 
CDD, the performance is worse. While \if) y ) is the satu- 
rating component, the fidelities for the other two states 
are lower. Notably, at the first concatenation level the 
fidelity is smaller than at the other levels. Relations sim- 
ilar to Eq. fl5B) hold for XZ CDD as 

^jXZPDDj^ ^ _ £jXZCDD2 _ ^yXZCDD3 _ _ ^jXZCDDll 



Furthermore, since in our case the in-plane axis errors, 
m x and n y , are negligible, pulse errors actually contribute 
only in the 3rd and higher orders. The good performance 
of SDD in Figure [3] compared to XY PDD is hence ex- 
pected. The SDD sequence we study here has a similar- 
ity in spirit with the Eulerian decoupling (EDD)£L. By 
making the structure of the decoupling sequence more 
symmetric, both SDD and EDD cancel the pulse errors 
to first order. 

Recall that in XY PDD, spin component \tp z ) is pre- 
served against pulse errors for a much longer time than 
the other two components. For SDD, such asymmetric 
behavior is less obvious. This is because to 3rd order, 
L/xy-sdd is a rotation operator with the rotation axis in the 
x-y plane (the expression is cumbersome and not shown 
here). Symmetrized DD treats the three spin components 
on a more equal footing. 

In the symmetrized sequence, Eq. ©, two 7TY-pulscs 
are adjacent to each other in the middle of a period. Since 
tjV/jV — i ; one would expect that removing these two im- 
perfect pulses from the sequence will yield better fidelity. 
However, the concatenated decoupling sequence based on 
such a reduced version of the SDD behaves much worse 
than the full sequence, as seen in Figure 0] This can also 
be understood from the error dependence of the evolu- 
tion operator. While [/ XY - SDD is identity up to first order 



10 




10 20 
N (number of cycles) 



FIG. 3: (color online). Simulation results for XY-based SDD. 
Fidelity as a function of the number of cycles is shown for 
different initial states. Fidelities in XY PDD are plotted as 
broken lines for comparison. Note that the number of pulses 
in one cycle of SDD is twice as large as that in PDD. 



in the pulse errors, with the two 7rY-pulscs removed, the 
evolution operator involves e y to first order 



^reduced 1 ct y <J 



(29) 



where the subscript "reduced" denotes that all adjacent 
7TY-pulses are removed from the sequence. Therefore, 
these two -ny-pulscs help keep a balanced structure for 
the symmetrized XY sequence to suppress error accumu- 
lation. 
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FIG. 4: (color online). Simulation results for concatenated 
XY-based SDD. Fidelity as a function of the concatenation 
level is shown for different spin states, (a) Results for the 
full sequence, (b) Results for the reduced sequence (with all 
adjacent 7TY-pulses removed). 



IV. CONCLUSIONS 

We have studied the effects of pulse imperfections on 
dynamical decoupling. Using the electron spins of phos- 
phorus donors in silicon as a testbed, we investigated the 
performance of different DD protocols. Experimental re- 
sults are qualitatively explained by theoretical analysis 
and quantitatively reproduced by numerical simulations. 

Dynamical decoupling sequences are susceptible to 
pulse imperfections. The accumulation of small pulse er- 
rors can have appreciable influence on the performance of 
DD. In two-axis based PDD, the spin component along 
a certain direction is preserved with high fidelity for a 
long time while the components in the perpendicular di- 
rections decay rapidly. DD sequences are analyzed in 
terms of spin rotation, and each DD sequence has an ef- 
fective axis in the presence of pulse errors. During the 
decoupling, viewed stroboscopically, the spins perform 
rotations about this axis. The better preservation of a 
particular spin component is attributed to the fact that 
this component is parallel to the effective axis of the DD 
sequence and hence is a quasi integral of motion, in spite 
of pulse errors. 

The XY-based decoupling protocols prominently out- 
perform their XZ-based counterparts. The different per- 
formance is explained analytically in terms of pulse er- 



rors. A symmetrized DD sequence suppresses the accu- 
mulation of first order pulse errors systematically and 
thus performs better than regular PDD. Interestingly, 
the accumulation of pulse errors is better suppressed by 
keeping the seemingly trivial adjacent identical pulses, 
7TY7TY7 in SDD. The symmetrized sequence itself has an 
error-balance structure in which these two pulses are in- 
dispensable. 

The concatenated dynamical decoupling sequences are 
found to be error-resistant. Although the number of 
pulses increases exponentially with the concatenation 
level, CDD have exactly the same error effects up to 
first order in the pulse errors. Concatenation can also 
alleviate the dependence of the error accumulation on 
the inter-pulse delay as well as the inhomogeneity in the 
static magnetic field. XY CDD, with the advantage of 
both the slow accumulation of pulse errors compared to 
XZ-based sequences and the superior error-resistance of 
the concatenation structure, demonstrates the ability to 
store all three spin components for a long time. 

In a general situation, when the internal bath dynam- 
ics is important, the joint action of the fluctuating bath 
and the imperfect pulses leads to a very complex qubit 
dynamics. However, in the case where the two effects 
are small, we can, at least on a qualitative level, consider 
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their contributions separately^ 4 -. As such, our work, fo- 
cused on the role of the pulse error accumulation, is com- 
plementary to the studies of the bath dynamics and the 
resulting homogeneous dephasing. By treating the bath 
field as static, we can single out the effect of the pulse 
error accumulation, and clarify its important role in the 
dynamical decoupling experiments. Moreover, the inho- 
mogencous static broadening is often important by it- 
self. For instance, in the experiments using single qubits 
with the projective readout (such as most experiments 
on quantum dots, the nitrogen-vacancy centers in dia- 
mond, superconducting qubits, etc.), where many exper- 
imental runs are needed in order to build statistics, the 
quasistatic shot-to-shot variations lead to inhomogeneous 
dephasing, and have to be refocused. 
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Appendix A: Radiation Damping 

In NMR and ESR experiments, when the in-plane mag- 
netization signal is strong, its precession about the static 
field induces a noticeable time-dependent current in the 
resonator. This current generates an ac magnetic field, 
which exerts a torque on the spins, and causes them to 
relax back to their equilibrium states. This phenomenon 
is called radiation damping (RD)i 42 i 43 We now consider 
the effect of RD on dynamical decoupling in our ESR 
experiments. 

The magnetic field generated by RD (the RD field), 
B r , is in the x-y plane, with its magnitude proportional 
to the in-plane magnetization, M, and its direction per- 
pendicular to 

7 e B r = —(-Myx + M x y) , (Al) 

where M. is the overall magnetization normalized by the 
equilibrium magnetization, M 

M = M/M , 



and the damping time constant 

r r = (A2) 

is determined by the quality factor, Q, the filling factor, 
•q, of the microwave-frequency resonator, and the equi- 
librium magnetization, Mo, which is determined by the 
temperature, T, and the doping density of 31 P. 

In our experiment the static field is Bq = 3500 Gauss, 
temperature T — 8 K, the filling factor and the quality 
factor of the resonator arc r\ = 0.5, Q — 500, respectively. 
The doping density of 31 P is very low, 5 x 10 14 cm -3 , 
corresponding to a damping time constant r r ~ 500 fjs, 
which is much longer than and hence the radiation 
damping effect is negligible in our DD experiments. 

However, for a higher doping density, of the order of 
10 15 " 17 (cm) -3 , RD could be noticeable. We calculated 
the influence of RD on the performance of the DD se- 
quence with r r — 2 fjs (corresponding to a doping den- 
sity ~ 1.5 x 10 17 cm~ 3 and other parameters the same 
as in the experiments). We find that in the dynamical 
decoupling process with imperfect pulses, RD could in- 
duce an asymmetry to the evolution of the initial states 
corresponding to the electron spin along the +z and — z 
directions. 

In Sec. IIII1 when we studied the effect of pulse errors 
on the DD performance, we neglected the finite pulse du- 
ration and treated the pulses as instantaneous rotations. 
RD, however, has its largest impact on spins during the 
pulses. If the pulses are ideal, the RD effects that take 
place during the delays before and after each pulse can- 
cel each others 3 - Moreover, here we focus on the initial 
states parallel or anti-parallel to z-axis. For a single elec- 
tron spin during the inter-pulse delays, the in-plane com- 
ponent of the spin magnetic moment is nonzero due to 
pulse errors (but is small). For a large number of elec- 
tron spins, these nonzero components are expected to be 
evenly distributed in the x-y plane, giving almost zero 
net in-plane magnetization. The magnitude of the RD 
field, B r , is proportional to the in-plane magnetization, 
hence little RD can be induced. During the pulse, on 
the other hand, the in-plane magnetization has nonzero 
value, and RD has finite contributions. 

We therefore take into consideration the finite duration 
of the magnetic field pulse. The pulse is implemented by 
applying a driving field to the electron spin for a finite 
duration, t p . For the spin state along the +z or — z direc- 
tions, at the edges of a pulse the in-plane magnetization 
is still small, and the resulting RD is negligible. This is in 
contrast with the pulse errors analyzed in Sec. IIII1 where 
contributions from the pulse edges were important. So 
here we can assume that the driving field for the pulse 
is turned on and off abruptly, and during the pulse the 
magnitude of the alternating field is constant for each 
spin. That is, in the rotating frame the driving field as 
a function of time is a square wave. In simulations in 
the rotating frame the driving field of the pulse 7Tx(7Ty) 
is taken to be along the axis n(m). The amplitude of the 
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driving field is taken as 

B p = B p + SB p (A3) 
where the average B p = 7r/(7 e i p ) and the deviation 

SB p = e/( le t p ) (A4) 

with e drawn from distribution Eq. (|13p . The parame- 
ters for pulse errors (angle errors e x , e y and the rotation 
axis errors n z , m z ) have the same magnitudes and dis- 
tributions as specified in Section III Dl During the pulse, 
both the static field and the driving field are present. So 
in the rotating frame the Hamiltonian, Eq. ([2]), also par- 
ticipates in driving the evolution during the pulse. We 
took the magnitude of the driving field B p = 1 Gauss 
as in experiment and correspondingly the duration of a 
7r- pulse is t p = 0.18 /xs. 

We simulated XY CDD with concatenation levels up 
to 4. for the electron spin initial states along the +z and 
—z directions. The fidelities are denoted as F +z and F_ z 
respectively. Here —z refers to the equilibrium state of 
the spin in the static field, and the +z state is prepared 
by applying a preparatory 7rx-pulse to the spin in the — z 
state. Due to the imperfection in the preparatory pulse, 
the fidelity for +z state becomes slightly different from 
the fidelity of the -z state. See Table U (A). 

The difference between states +z and — z produced 
by the preparatory pulse is further amplified by the RD 
in the concatenated DD sequence. Figure [5] shows the 
simulation results for XY CDD. The decoupling protocol 
preserves the state —z with higher fidelity than +z. The 
difference is more prominent for concatenation level 4. 

This is because RD only takes effect during the mag- 
netic field pulses when the in-plane magnetization is 
nonzero. The RD time constant, r r , is larger than the 
pulse duration by two orders of magnitude. The effect of 
RD is therefore very small for the first three concatena- 
tion levels, and after the accumulation during hundreds 
of pulses, RD shows up prominently at level 4, (FigureEJ). 

To better understand the different roles of the prepara- 
tory pulse, the effect of RD during the pulses and dur- 
ing the inter-pulse delays, simulations are performed for 
three cases: (A) without RD, i.e. r r — oo; (B) RD is 
only included during inter-pulse delays; (C) RD is in- 
cluded throughout the whole evolution (i.e., both during 
and between the pulses). Table Q] shows the data for XY 
CDD in the three cases. Compared to the results with- 
out RD, the difference between F +z and F- z for levels 



1-3 is mainly due to the preparatory pulse, and is little 
affected by RD. For level-4, RD contributes appreciably. 
The different fidelities for level-4 in cases (B) and (C) 
confirms that RD contributes mostly during the pulses. 

Conventionally, starting from the linear Bloch equa- 
tions, the states +z and — z would have exactly the 
same dynamical behavior. Radiation damping breaks 
this symmetry. The RD field, B r (Eq. (|Al|) h is a macro- 
scopic field and is dependent on the overall magnetization 
of the electron spin ensemble. The Bloch equation is then 




1 1 1 1 1 

12 3 4 

Concatenation Level 



FIG. 5: Fidelities F +z and F- z for XY CDD sequence with 
radiation damping considered. B p = 1 Gauss and r r = 2 /is. 



TABLE I: XY CDD pulse sequence. Dependence of the fideli- 
ties on the concatenation level for electron spin initial states 
along the +z and —z directions. Pulse field B v — 1 Gauss. 
Results for three cases are shown: (A) r r = oo; (B) r T — 2 /is, 
RD takes effect only during inter-pulse delays. (C) r r = 2 fis, 
RD takes effect all the time. 



level 


(A) F +z 


F- z 


(B) F +z 


F- z 


(C) F +z 


F- z 


1 


0.929 


0.999 


0.930 


0.999 


0.929 


0.999 


2 


0.923 


0.994 


0.923 


0.995 


0.924 


0.993 


3 


0.928 


0.998 


0.926 


0.998 


0.920 


0.994 


4 


0.924 


0.994 


0.877 


0.933 


0.668 


0.832 



no longer linear in the magnetization, M, resulting in the 
different dynamics of the quantum states along the +z 
and — z directions. 
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